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Using spin-spiral technique within the full potential linearized augmented-plane-waves (LAPW) 
electronic structure method we investigate the magnon spectrum and Neel temperature of bcc 
Eu. Ground state corresponding to an incommensurate spin-spiral is obtained in agreement with 
experiment and previous calculations. We demonstrate that the magnetic coupling is primarily 
through the intra-atomic / — s and f — d exchange and Ruderman-Kittel-Kasuya-Yosida mechanism. 
We show that the existence of this spin-spiral is closely connected to a nesting feature of the Fermi 
surface which was not noticed before. 

PACS numbers: 71.18.+y, 71.27. +a,75.10.Hk 



I. INTRODUCTION 

The magnetic behavior of most rare-earth (RE) mate- 
rials is governed by localized magnetic moments interact- 
ing indirectly through the sea of delocalized valence elec- 
trons. A typical feature of RE systems is existence of two 
electron species, localized 4/'s exhibiting atomic like be- 
havior with strong Coulomb interaction and delocalized 
5d and 6s electrons with interaction merely renormaliz- 
ing the band dispersion. Yet, a typically weak interac- 
tion between these two species, either due to intra-atomic 
exchange or band mixing (hybridization), gives rise to 
a variety of magnetic behaviorsi. On the model level 
this behavior is captured by periodic Anderson model 
or Kondo lattice model^. The ab initio electronic struc- 
ture methods based on density functional theory (DFT) 3 
and the standard semi-local approximations^ have no- 
torious problems in dealing with the strong correlations 
within the 4/ shell. In particular the splitting between 
occupied and unoccupied 4/ bands, which is the way a 
single-particle band structure can capture the effect of 
Coulomb interaction, is missing. Consequently, the / 
bands appear at the Fermi level resulting in unrealis- 
tic filling of both / and valence orbitals. Two reme- 
dies can be used: (i) open-core treatment or (ii) addi- 
tional Coulomb term with the simplest example being 
the LDA+U method. In the open-core treatment the 4/ 
orbitals are kept separate from the rest of the valence 
Hamiltonian and the interaction with the valence states 
is only through the self-consistent potential. Obviously 
the proper filling of both / and valence bands is easy to 
achieve if integer, however all kinematic exchange effects 
(e.g. superexchange) based on band mixing are miss- 
ing. In the LDA+U approach the splitting between the 
occupied and unoccupied / band is obtained due to an 
additional orbital-dependent term. All possible exchange 
processes are in principle accounted for in this approach. 

Compounds containing Eu in 2+ formal valency are 



particularly well suited for LDA+U treatment, since the 
orbital degrees of freedom are quenched in the half filled 
/ shell. Examples involve ferromagnetic insulators EuO 
and EuS&, semimetal EuBgi and metallic elemental Eu. 
Presumably the exchange mechanisms in these materials 
are quite different. 

In this paper we investigate the magnetic ground state 
and magnon spectrum of elemental Eu. Europium crys- 
tallizes in body-centered cubic (bcc) structure with a 
lattice constant^ of 4.555 A at 100 K. The magnetic 
groundstate was found to be a spin spiral and the Neel 
temperature^ of 91 K. The electronic structure was pre- 
viously investigated by Freeman and DimmockiS and 
Andersen and Loucksii using the Xa potential. Re- 
cently Turek et alH used a real-space perturbation ap- 
proach based on tight-binding linear muffin-tin orbital 
(TB-LMTO) method to calculate the exchange parame- 
ters and corresponding magnon spectrum and Neel tem- 
perature. Here we use a reciprocal space based spin spiral 
approach which can be viewed as complementary to the 
real space calculations. Unlike the above authors who 
used the open core treatment of the 4/ orbitals we employ 
the LDA+U method. We use the linearized augmented- 
plane-waves (LAPW)i£ method and its extension to non- 
collincar magnetic structures utilizing the generalized 
Bloch theorem^ for calculation of the spin spiral states. 
We interpret our results in terms of Ruderman-Kittel- 
Kasuya-Yosida (RKKY) exchange mechanismi£ii£ii£ and 
the Fermi surface property. For this purpose we evaluate 
the low frequency limit of the generalized susceptibility 
and identify the nesting features of the Fermi surface. 



II. COMPUTATIONAL METHOD 

We have used the Wien2ki& implementation of the 
full potential LAPW method and its extension for non- 
collinear spin structures^. The effective single-particle 



FIG. 1: (color online) The band structure as obtained with 
LDA+U method. The brighter (red) lines correspond to ma- 
jority spin. In the side panels the corresponding densities of 
state for the minority (left) and majority (right) spin projec- 
tions are shown. 



FIG. 2: The q-dependent exchange parameter calculated us- 
ing the spin-spiral approach. Comparison of the results for 
different k-point samplings indicates that the 15x15x15 mesh 
is reasonably well converged. 



III. RESULTS AND DISCUSSION 



potential with constructed from LDA+U functional, with 
the exchange-correlation potential of Perdew and Wang 5 
and the double-counting scheme of Anisimov et alm^. 
The Coulomb term parametrized with U (7 eV unless 
stated otherwise) and J (0.75 eV) was applied to the 4/ 
orbitals. The spin-spirals were treated using the gener- 
alized Bloch theoremiA which prohibits inclusion of the 
spin-orbit coupling. The atomic sphere approximation 
was employed, in which the direction of the exchange 
field is constrained inside atomic spheres and allowed to 
vary continuously in the interstitial space. 

A spin spiral is defined by a propagation vector q and 
angle 9 between the local moment and the precession 
axis. The orientation of the precession axis itself is arbi- 
trary unless the spin-orbit coupling is taken into account. 
Each spin spiral was calculated selfconsistently. This ap- 
proach has two deficiencies compared the perturbative 
approach employing the force theorem£i:(i) it costs much 
more computational effort, (ii) more importantly one has 
to work with total energies instead of the sums of eigen- 
values, i.e. looking at small differences of large numbers 
which requires high accuracy. As for (i), when perform- 
ing calculations for q-vectors along a certain path in the 
reciprocal space converged spin density from a nearby q 
can be used as starting point, which reduces the num- 
ber of iterations needs significantly compared to starting 
from scratch for each q. As for (ii) , dependence of spin 
spiral energies can be studied without being limited to 
small angles. 



A. Bandstructure 

In Fig. ^ we show the spin-projected bandstructure 
obtained with U of 7 eV and J of 0.75 eV. The lowest 
valence band around the T point with a predominant 6s 
character corresponds to a strongly dispersive s-band. 
Moving toward the zone boundary mixing with e?-band 
takes place and the doublets at H and P points have a 
pure d symmetry. The states in the vicinity of the Fermi 
level have mostly d character. 

The occupied 4/ levels are localized about 3 eV below 
the Fermi level and cross the lowest valence band, with 
negligible hybridization, which is reflected by completely 
flat dispersion. The unoccupied 4/ bands are approx- 
imately 7 eV above the Fermi level. A bandwidth of 
about 2 eV originates from mixing with the 6p and 5<i 
bands. Lack of mixing with valence states in the oc- 
cupied / bands indicates that the kinematic exchange, 
involving hopping from localized / orbitals into the de- 
localized band states, is not important here. The inter- 
action between the localized / states and the rest of the 
electronic system, is dominated by intra-atomic / — s and 
/ — d exchange. Consequently the spin polarized bands 
and density of states below 2 eV exhibit almost a perfect 
rigid shift. The deviations from this pattern are mostly 
due to the splitting of the s-like bands being less than 
that of d-like bands, because the more extended s orbital 
has smaller overlap with the polarized / density. 
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FIG. 3: The q-dependent exchange parameter renormalized 
to satisfy the sum rule J5J along the high symmetry lines; ab 
initio data (symbols), Fourier interpolation (line). The right 
panel shows the corresponding density of states. 

B. Magnon spectrum 

The magnetic excitations are discussed in terms of a 
classical Heisenberg Hamiltonian 

H = - ^2 JRR/e R • e R /. (1) 

RR' 

A spin spiral characterized by the propagation vector q 
and angle 9 has the form 

e R = (sm(0) cos(q • R),sin(0) sin(q • R). cos(0)), (2) 

The corresponding energy per lattice site obtained from 

E(q) = sin 2 (0)(j(q) - J(0)) + J(0) + E (3) 
J(q) = Joi? cxp(iq • R). (4) 

R 

is to be compared to the ab initio results. Eq is the non- 
magnetic part of the total energy. Note that only the 
difference J(q) — J(0) can be obtained from the knowl- 
edge of E(q). In order to fix the value of J(q) the sum 
rule 

J dqj(q) = 0, (5) 

where the integration is over the Brillouin zone, is to be 
employed. The J(q) normalized this way contains only 
information about the inter-site Jrr' and can be related 
to the ordering temperature, which we discuss below. 

Using the Hamiltonian ^ involves several approxima- 
tion. Quenching of the orbital moment in a half filled 



FIG. 4: (color online) The paramagnetic ab initio bands 
(dots) together with the smooth Fourier interpolation of the 
low valence bands (red lines). Note that the inaccuracies of 
the fit at the band crossings (for example midway between P 
and H) are far from the Fermi level. 

shell, rigidity of the / moment and its size (S=7/2) well 
justify the use of Heisenberg Hamiltonian in classical ap- 
proximation. In addition we assume that the exchange 
parameters are constant. This is not a priori guaran- 
teed, since the electronic structure of the band electrons, 
which carry the exchange interaction between the local / 
moments, depends on the arrangement of local moments. 
This question was discussed in detail by Nolting et al22k 
who derived an expression for the effective exchange pa- 
rameters Jrr' in terms of the conduction electron self- 
energy. Experimentally this leads to temperature depen- 
dence of the effective exchange parameters. If this effect 
were important a deviation from the 9 dependence of eq. 
© is expected. The fact that we have not found any sig- 
nificant deviation from (|3J) in the range from 90° to 30° 
serves as a justification of use of Hamiltonian (T|. This is 
also agreement with a rather big ratio of the band- width 
to the exchange splitting in the conduction band. 

In Fig. |2 we show the q-dependent exchange parame- 
ter obtained with maximum 9 of 90°. Calculations per- 
formed with U of 6 eV and 8 eV lead to almost identical 
dispersion supporting our previous conclusion about the 
f — d exchange mechanism which does not depend on the 
position of the / bands. The calculations yield a mini- 
mum energy corresponding to a spiral with propagation 
vector Q of about 0.3 x (^,0,0) and another two local 
minima along TN and TP lines. ' In order to address the 
Neel temperature the E(q) throughout the Brillouin zone 
is needed. To this end we have calculated the spin-spirals 
on a 10 x 10 x 10 regular q-grid (44 irreducible q-points) 
and used a smooth Fourier interpolation 2 ^^ 2 ^ to obtain 
J(q) on a denser grid. In Fig. |3| we compare the inter- 
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FIG. 5: The paramagnetic Fermi surface: the lobed tetrahe- 
dra are centered at 8 P points and 'supereggs' at 6 H points on 
the surface of the 1st Brillouin zone (only one representative 
of both species is shown for sake of clarity) . 



polating function to the ab initio results along the high 
symmetry directions. Since only a few high symmetry q- 
points were contained in the regular grid the agreement 
between the interpolating function and the ab initio data 
points indicates the quality of the grid. Comparison to 
the results obtained by Turek et at reveals a very good 
agreement of the key features. Besides the similar over- 
all shape of the spectrum, we have obtained the same 
ordering of the peak according to their size as well as 
similar absolute values. The most apparent deviation 
is the position of the ferromagnetic state J(0) which is 
relatively more favorable in our calculation. Evaluating 
the expressions^ for Neel temperature in the mean-field 
(rpMF^ anc j ranc i om _ph ase (T$ PA ) approximations 




FIG. 6: (color online) The nesting function in the FHNH 
plane. The lower-left and upper-right corners correspond to 
r point with the typical 1/q divergence, the remaining corners 
are at H point and the N point is in the center of the square. 
The upper panel shows the sum over all (2) bands. The mid- 
dle panel shows the 'superegg' contribution and the lower 
panel the lobed tetrahedra contribution (the contribution of 
the superegg-to-lobed tetrahedron processes is not shown). 



k B T* 1F = |j(Q) (6) 

{k B Tr A r x = ~ e{ t j (Q) - j (q)r l + 

[J(Q)-^(q + Q)-^(q-Q)] -1 } 

(7) 

with the interpolating function we arrive at T^ F =170 K 
and Tft PA =U2 K which again agree very well with 147 
K and 110 K obtained by Turek et at. as well as with 
the experimental value 8,9 of 90.5 K. 



C. Fermi surface and nesting 

In this section we discuss the paramagnetic Fermi sur- 
face and its connection to the spin-spiral ground state. 
Very fine k-point sampling is required in order to study 
the fine details of the Fermi surface. We adopted the fol- 
lowing approach. Starting with 100 irreducible k-points 
obtained with LAPW code we have used the procedure 
for smooth Fourier interpolation^ for the bands at the 
Fermi energy. We have generated the band energies on 
a finer mesh of approximately 1.5 x 10 6 k-points in the 
whole Brillouin zone, which was used for the calculations 
reported below. In Fig. ^ we show the paramagnetic 
bandstructure together with the Fourier interpolation, 
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which is excellent near Ep (note that the interpolation 
was performed on different-regular mesh). 

The paramagnetic Fermi surface of bcc Eu, Fig. 03 
consists of an electron pocket centered at H point, and a 
hole pocket located at P point. Symmetry related degen- 
eracy of the valence band along the P-H direction results 
in touching of these pockets. The first numerical investi- 
gation of the Fermi surface of bcc Eu was reported by An- 
dersen and Loucksii, who called these pockets 'superegg' 
and 'tetracube' respectively. While our calculation repro- 
duces the 'superegg' shape, contrary to Andersen's cube 
with lobes we find rather a rounded tetrahedron with 
lobes. 

In order to make connection between the Fermi surface 
geometry and calculated spin spiral dispersion we have 
evaluated the imaginary part of the generalized suscep- 
tibility 



Xij(q>w) = - 



/( ei (k))-/(e,(k + q)) 



(8) 



where e(k) is the band energy, /(e) is the Fermi-Dirac 
function and i and j are band indices. In the limit lo — > 
the imaginary part of x(q, cj) behaves as 



Imx y (q,w) 



7iWy(q), 



e F )S{e 3 (k + q) - e F ) 



(9) 
(10) 



where v(c\) measures so called nesting, i.e. the extent 
to which different parts of the Fermi surface weighted 
by the inverse square of the Fermi velocity are parallel. 
While the real part of the generalized susceptibility is 
directly related to the exchange parameters J(q) it is the 
imaginary part which has a straightforward geometrical 
interpretation. The real and imaginary parts are bound 
by the Kramers-Kronig relations, which for lo — limit 
read 



Re X (q,^ = 0) = i^ 



(11) 



Comparing this to one can readily see that a 
peak in ^(q) provides a significant contribution to 
a peak in Rex(q, lj = 0). In the special case 
q — > limit z/(q) diverges as l/\q\, however the limit 
lim 9 ^ lim^^o R e x(q: w ) remains finite, equal to the 
density of states (or corresponding partial density of 
states) at the Fermi level N(cf)- In general peaks in 
i^(q) indicate a tendency toward instability of the Fermi 
surface toward formation of incommensurate structures 
such as spin or charge density waves. 

Andersen and Loucks concluded that the origin of the 
spin spiral groundstate is the nesting between the oppo- 
site faces of their 'tetracube'. However, our investigation 
provides a different picture. In Fig. [|J]we show the nesting 
function ^(q) in the rHNH plane. There is a prominent 




FIG. 7: (color online) The nesting function in the FHPN 
plane. The arc shaped features intersecting the TP line orig- 
inate from superegg-to-lobed tetrahedron processes. 



peak on the T-H line and a weaker feature on the T-N 
line. In order to identify the origin of these features we 
have calculated separately the contributions of the differ- 
ent sheets of the Fermi surface. Fig. [^demonstrates that 
both peaks originate from transitions between the lobed 
tetrahedra surfaces. The q parallel to the r-H direction 
amounts to moving the lobed tetrahedron centered at the 
P point toward the next P point on the same face. The 
peak on the r-H line correspond to the overlap of the 
lobes, which is illustrated explicitly in Fig. [S] The peak 
on the T-N line corresponds to the sum of two nesting 
vectors along adjacent T-H directions. Therefore it orig- 
inates from an overlap of the tips of the lobes, however, 
belonging to tetrahedra centered at P points which are 
not on the same face. In Fig. we show the nesting 
function in the THPN plane. Besides the peaks on T-H 
and T-N lines we find weaker features on the T-P line. 
Analysis of the contributions from different sheets shows 
that these features originate from transitions between the 
two types of Fermi surfaces. 

Eventually we can compare the information obtained 
from the nesting function to the magnon dispersion. We 
have found a prominent nesting feature on the T-H line 
and weaker one on the T-N line. These are related to the 
minima in the spin-spiral energies. The energy minima 
do not sit exactly at the positions of corresponding the 
nesting vectors. While the nesting vectors contain in- 
formation about the states directly at the Fermi energy, 
all states contribute to the spin-spiral minima, yet with a 
weight decreasing with the distance from the Fermi level. 
A discussion of the generalized susceptibility in the con- 
text of electron-phonon coupling can be found in Ref. 
|2^.The features on the T-P line are rather weak and do 
not fit well with the position of the T-P peak therefore 
we do not draw any conclusion about its direct relation 
to the local energy minimum on the T-P line. 
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FIG. 8: The lobed tetrahedron translated by q = 0.4 x 
(2^,0,0) corresponding to the F-H peak elucidates the ori- 
gin of the spin-spiral groundstate. 

IV. CONCLUSIONS 

Using a full-potential method with LDA+U functional 
and spin-spiral approach we have obtained a magnon 



spectrum which is in good agreement with that published 
by Turek et al. using TB-LMTO. The calculations repro- 
duce well the experimentally observed spin-spiral ground- 
state and provide a reasonable estimate of the Neel tem- 
perature. The good agreement of two rather different 
computational methods indicates a robustness of these 
physical properties. Moreover we have shown that the 
values of the exchange parameters are insensitive to the 
precession angle of the spin spiral as well as the value 
of U, i.e. exact position of the occupied / bands. This 
confirms the picture of bcc Eu as a Kondo-lattice sys- 
tem with ferromagnetic exchange of intra-atomic origin 
in the RKKY regime. We have identified the origin of the 
spin-spiral ground state in terms of nesting properties of 
the Fermi surface. In particular we have shown that the 
nesting is connected to the lobed tetrahedron hole pock- 
ets centered at the P point of the Brillouin zone. 
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